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Abstract 

A competitive growth model (CGM) describes aggregation of a single type of particle under two 
distinct growth rules with occurrence probabilities p and 1 — p. We explain the origin of scaling 
behaviors of the resulting surface roughness with respect to p for two CGMs which describe random 
deposition (RD) competing with ballistic deposition (BD) and RD competing with the Edward 
Wilkinson (EW) growth rule. Exact scaling exponents are derived and are in agreement with 
previously conjectured values Q, Q]. Using this analytical result we are able to derive theoretically 
the scaling behaviors of the coefficients of the continuous equations that describe their universality 
classes. We also suggest that, in some CGM, the p— dependence on the coefficients of the continuous 
equation that represents the universality class can be non trivial. In some cases the process cannot 
be represented by a unique universality class. In order to show this we introduce a CGM describing 
RD competing with a constrained EW (CEW) model. This CGM show a transition in the scaling 
exponents from RD to a Kardar-Barisi-Zhang behavior when p — > and to a Edward Wilkinson 
one when p — ► 1. Our simulation results are in excellent agreement with the analytic predictions. 

PACS numbers: 81.15.Aa, 89.75.Da, 68.35.Ct, 05.10.Gg 
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I. INTRODUCTION 



Evolving interfaces or surfaces are of great interest due to their potential technological 
applications. These interfaces can be found in many physical, chemical and biological pro- 
cesses. Examples include film growth either by vapour deposition or chemical deposition P| , 
bacterial colony growth 4] and propagation of forest fire (f| . 

For a system exhibiting dynamical scaling, the r.m.s. roughness W of an interface is 
characterized by the following scaling with respect to time t and the lateral system width 
L: 

W(L,t)~L a f(t/L z ), 

where the scaling function f(u) behaves as f(u) ~ viP with (3 = z/a for « C 1 and f(u) ~ 
constant for u> 1. The exponent a is the roughness exponent that describes the dependence 
of the saturated surface roughness versus the lateral system size, while the exponent (3 
describes scaling at an early stage when finite-size effects are negligible. The crossover time 
between the two regimes is t s = L z . 

A widely studied phenomenological equation representing the non-equilibrium growth of 
such interfaces is the Kardar-Parisi-Zhang (KPZ) equation. In 1 + 1 dimensions, it states 
that: 

dhix.t) d 2 h(x,t) , ( dh(x,t)\ , s 

where h(x, t) is the local surface height at lateral coordinate x and time t. The coefficients v 
and A represent the strength of the linear and nonlinear surface smoothing terms respectively. 
The noise f]{x,t) is Gaussian with zero variance and covariance 

(r](x, t) v(x',t')) = 2D 5(x - x') 8(t - t') 

where D is the strength of the noise. The exponents characterizing the KPZ equation in 
the highly nonlinear strong coupling are a = 1/2 and /? = 1/3. In contrast, at A = the 
linear Edward Wilkinson (EW) equation is recovered leading to the weak coupling exponents 
a = 1/2 and (3 = 1/4. When both v and A are zero the growth reduces to simple random 
deposition (RD) with (3 = 1/2 but a lack of any saturation regime. 

There has been a recent interest in the study of competitive growth models (CGM) 
analyzing the interplay and competition between two current growth processes for a single 
surface. These CGM are often more realistic in describing growing in real materials |5|, 



in which more than one microscopic growth mode usually exist. For example, two distinct 

Erowth phases were observed in experiments on interfacial roughening in Hele-Shaw flows 
3,0] as well as in simulations on electrophoretic deposition of polymer chains B, |icj|. The 
resulting universalities from these competing processes are not well understood (jj, |2|, [ll|, [12 1 . 
Recently Horowitz et al. Q] introduced a CGM called BD/RD in which the microscopic 
growth rule follows either that of the ballistic deposition (BD) model with probability p or 
simple random deposition (RD) with probability 1 — p. This system exhibits a transition at 
a characteristic time from RD to KPZ. They found numerically that the scaling behavior of 
W in that model is giving by the empirical form 

W ~ —jF ( — ( 2 ) 
p s \ v -vL z ) y ' 

Based on numerical estimates, exact values 

5 = 1/2 and y = 1 (3) 

have been conjectured for the BD/RD model. Based on this conjecture, the authors con- 
cluded using scaling arguments that the model follows Eq. (0) with v ~ p and A ~ p 3 ^ 2 . 
Subsequently, a similar CGM namely EW/RD describing the competition between the EW 
model and RD Q. The simulations showed that Eq. © also holds with a different set of 
exponents which are conjectured to be 

5 = 1 and y = 2 (4) 

This model can also be described by Eq. (0) with v ~ p 2 and A = ll. 

In this paper, we explain the scaling form (J2J) and derive rigorously the exact exponents S 
and y using simple arguments. In addition, from the above examples of CGM, one might be 
tempted to conclude that a CGM based on RD and a model in the KPZ (EW) universality 
class should always lead to an overall process in the KPZ (EW) class. We suggest that these 
naive predictions of the universality is not always correct. Close examination of the micro- 
scopic details of the growth models is indeed essentially. This is illustrated by introducing a 
constrained EW (CEW) model. Although this model essentially belongs to the EW class, a 
CGM in the form CEW/RD at sufficiently small p results in an overall process in the KPZ 
universality class. It also demonstrates that a CGM can crossover from one universality 
class to another by varying p. 



II. EXACT SCALINGS FOR CGMS 



In the RD model a particle is dropped at a randomly selected column increasing the 
local surface height by one. For the CGMs BD/RD or EW/RD described above, a particle 
is deposited on the surface following a RD process with probability 1 — p and by another 
process A (which is either BD or EW) with probability p. Now, we derived analytically the 
exact exponents 5 and y which characterize the p dependence of the scaling behavior of W 
given in Eq. (J2J). We consider p — > 0. At each unit time, L particles are deposited. The 
average time interval between any two consecutive A events at any column i is r = 1/p. 
During this period, r — 1 ~ r atoms on average are directly stacked onto the surface 
according to the simple RD rule. The local height at column i hence increases by r\i which 
is an independent Gaussian variable with mean Tj = r and standard deviation a v = y/r 
according to the central limit theorem. The mean however only leads to an irrelevant rigid 
shift of the whole surface. We can easily apply a vertical translation so that rj = 0. After 
these t RD events at column i, one A event on average is expected at the same column. 

Now we consider A to be the BD process. The CGM is then the BD/RD model. When 
a BD event occurs at column i, its height is updated in the simulations according to hi — > 
max{hi-i, hi + i, hi + 1}. In the limit p — ■> so that a v >> 1, the height of the atom is 
negligible compared with the increments due to the RD events. The growth rule hence 
reduces to 

hi — > max{/ij_i, h i+ i, hi} (5) 

We have now arrived at a limiting BD /RD model defined as follow: At every coarsened time 
step r = 1/p, the local height hi at every column i first changes by an additive Gaussian 
noise term rji with mean zero and standard deviation a v = \fr. Then the limiting BD growth 
rule in Eq. (jSJ) is applied to every column i. A more careful analysis should account for 
the fact that the BD events at various columns indeed occur randomly and asynchronously 
during the period r but this will not affect our result. In this limiting model, the time and 
the vertical length scales are determined completely by r and a v respectively. Therefore, 
time scales ast~r = 1/p while roughness scales as W ~ a v ~ 1/p 1 ' 2 . This explains the 
p dependence of the scaling form in Eq. (j2J). In particular, we obtain the exact exponents 
y = 1 and 5 = 1/2 in agreement with values in Eq. (jSJ) first conjectured in Ref. |l( but not 
derived analytically before. 
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Next, we assume that A represents the EW growth rule instead. For this growth rule, 
a particle is dropped at a random column but when it reaches the surface it is allowed 
to relax to the lower of the nearest neighboring columns. If the heights at both nearest 
neighbors are lower than the selected one the relaxation is directed to either of them with 
equal probability. Our CGM now becomes a EW/RD model. The corresponding derivation 
of the characteristic length and time scales is similar to that for the BD/RD case. Assuming 
again p — > 0, the average time interval between any two consecutive EW events at any given 
site is r = 1/p. Consider a characteristic time nr. On average nr RD events occur at any 
given site leading height increments with a standard deviation y/nr. However, only n EW 
events take place. The resulting smoothing dynamics is such that a big step for instance 
will typically decrease in height by n. For scaling to hold, the two length scales y/nr and n 
have to be identical and we obtain r ~ n. The characteristic time scale considered is hence 
nr ~ 1/p 2 , while a characteristic length scale for the surface height is ^frvr ~ 1/p. We have 
hence derived y = 2 and 6 = 1 previously conjectured in Ref. 

III. CEW AND CEW/RD MODELS 

We now introduce the constrained EW model (CEW) which is a generalization of the 
EW model. In 1 + 1 dimensions, particles are aggregated by the following rules. We choose 
a site i at random among the L possible sites. The surface height hi at the selected column 
is increased by one if this height is lower than the values hi±i at the neighboring columns. 
Otherwise either h^\ or h i+ i, whichever smaller, is increased to 

hi±i = max{/i i± i + l,hi- c} (6) 

If hi _i = hi+i, either one will be updated with equal probability. Growth at i ± 1 physically 
represents the rollover of a newly dropped particle to a lower site nearby under the influence 
of gravity for example. In the original EW model, there is no limit in the vertical distance 
transversed during the rollover and the particle can in principle slides down a very deep cliff 
if one exists next to column %. This is unphysical if there is a finite chance for the sideway 
sticking of the particle to the cliff. In the CEW growth rule defined in Eq. (jUJ) , this vertical 
drop during rollover is limited to c by the process of sideway sticking. As c — > oo, it is 
easy to see that CEW reduces to the standard EW model. At c = 0, sideway sticking of 
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particles occurs frequently and CEW behaves similarly to BD, although the precise growth 
rules are different. In the rest of this paper, we put c = 4. As will be demonstrated by the 





FIG. 1: Scaling plot of W/L a as function of t/L z for c = 4 and different L values: L = 256 (o), 
L = 512 (□), L = 1024 (o). (a) For p = 0.02 we use the scaling exponents characteristic of the 
KPZ equation a = 1/2 and z = 3/2. The dashed line with slope 1/2 and the dotted-dashed line 
with slope 1/3 are used has a guide to show the RD and KPZ regimes respectively (b) Show the 
same as (a) but for p = 0.64, with a = 1/2 and z = 2, characteristic of the EW. behavior. The 
dashed line with slope 1/2 and the dotted-dashed line with slope 1/4 are used has a guide to show 
the RD and the EW regimes respectively. 



simulations presented later, the CEW model at c = 4 at practical length and time scales 
belongs to the EW universality class. Finally, we can define the CEW/RD model which is 
a CGM based on RD and CEW. Similar to the definitions of other CGM defined above, at 
each simulation step, the CEW growth rule is applied with probability p while a RD event 
occurs with probability 1 — p. Time t is then increased by 1/L. 

Now we present the simulation results. In Fig. G^a) an d Fig. E^b) we show the Log-Log 
plot of W/L a as function of t/L z for two limiting values of p. For p — > the behavior is 
consistent with the KPZ universality class with a = 1/2 and z = 3/2, while for p — > 1 the 
system behaves as predicted by the EW equation with a = 1/2 and z — 2. The initial 
regime corresponding to the RD deposition, with (3 = 1/2 does not scale with the system 
size. In order to show that the universality class depends on p, we compute f3 as function of 
t using successive slopes defined in Ref jl3l |. 

Figure 121 plots (3 as a function of log t for different values of p. At the beginning (3=1/2 
as expected for the initial RD regime. After this early regime, the system evolves either to 
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FIG. 2: Plot of the dynamic exponent (3 as function of log t for L = 8192 and different p values: 
p = 0.02 (o), p = 0.04 (□), p = 0.08 (o), p = 0.16 (A), p = 0.32 (+), p = 0.64 (*), and p = 1.0 
(x) showing the change in the behavior of (5 with time. The (3 values where computed over 100 
realizations. The arrows are used as guides to show the asymptotic exponents 

the KPZ class with (3 = 1/3 for p — > or to the EW class for p — > 1 with (3 = 1/4. For 
intermediate p values (after the RD regime), the system behaves as in the weak coupling 
of Eq. (^). It is easy to observe a transition from an EW to a KPZ for p > 0.32 while for 
p — > the system always belongs to the KPZ universality class. 

IV. GENERALIZED CONTINUOUS EQUATIONS AND SCALING 

As we show above the CEW/RD model has a transition from a KPZ to a EW as the 
tuning parameter p goes from to 1. Thus the nonlinear coefficient A of the KPZ equation 
has to vanish as p — > 1. In order to understand the functional form of X(p) we perform a 
finite size scaling analysis of the growth velocity based on |14j : 



where Av(L,t) = v(L = 1024, t) - v(L = 10, t) and v(L,t) = (dh/dt). The Av correction 
should go to zero when the nonlinear term A vanishes. Thus, using Eq. ((7|) we can determine 
how A in the KPZ equation depends on p. In Fig Olwe plot Av as function of p for fixed L 
to show the p dependence of A. From the plot we can see that the functional form for the 



Av(L) ~ XL 



—a 



for t > t 



s 



(7) 
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FIG. 3: Linear plot of A as function of p. The averages where taken over typically 500 realizations. 
The errors bars where computed using 5 independent sets of 500 realizations each. 



CEW/RD CGM is totally different from the scaling form A(p) ~ p 3 ^ 2 of the BD/RD model 
P). For small p values A has a power law dependence with p, while for p — > 1, X(p) — > 
with a fast decay. Now we proceed to generalize Eq. (JJJ for the CGM model applying the 
following transformation : h = h f(p) b a , x = b x and t = g(p) t b z . where b is the length 
transformation. As the interface evolution of these CGM are independent of b they can be 
described, after applying the transformations defined above, by the generalized continuous 
equation: 

^ =V(P) ^ +X(p) (^y +Dip)n{x , t) (8) 

with 

"(p) = "og(p) (9) 
Hp) = A f(p) g(p) (io) 

m - j$ (id 

where g(p) is related to the characteristic time scale when the correlations begins to dominate 
the dynamic of the interface and f(p) is related to the saturation length scale. Notice that 
for the EW universality class Ao = Using the exact results from Sect HH 

p 1 / 2 , for BD/RD and CEW/RD when p -> 0, 

f{P) ~ { ; (12) 

p, for CEW/RD when p — > 1, 
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FIG. 4: Log- Log plot of W f(p) as function of t g(p) for L = 8192 and different p values. In (a) 
p = 0.02 (o), p = 0.04 (□) and p = 0.08 (o), g(p) = p and f(p) = p 1 / 2 . In (b) p = 0.32 (o), p = 0.64 
(□) and p = 1.0 (o), g(p) = p 2 and f(p) = p. Notice the departure from the EW scaling behavior 
for p = 0.32 



and 



. p, for BD/RD and CEW/RD when p -»• 0, 

0(P) ~ ^ „ , ( 13 ) 

p 2 , for CEW/RD when p -»• 1, 



Thus Eqs. (^. (fTUjl and (fTT|) can be replaced by: 
i/(p) ~ < 



i/ P, for BD/RD and CEW/RD when p -> 0, 
z/ p 2 , for CEW/RD when p — ► 1, 



(14) 



. A p 3/2 , for BD /RD when p -+ 0, 
A(p) ~ ^ (15) 
0, for CEW/RD when p -»• 1, 

and D(p) ~ D for BD/RD and CEW/RD independant of p Thus, for the CGM belonging 
to the KPZ (EW) universality class the evolution of the interface is given by: 

dh_{v p + Ac p 3/2 (^^) 2 + Vi(t), for BD/RD and CEW/RD when p - 

d * I i/ P 2 + for CEW/RD when p -+ 1, 



The scaling behavior of VT[11] is then: 



, F ( g (p)X ox /Mt) for KPZ, 
Wf(p)/L a ~ { ^ W V 1/0 ^ v ' ' (17) 

F Up) i/o , for EW, 



9 



that after replacing f(p) and g(p) given by Eq. (|T2*|) and Eq. (JT5)l leads to the exact scaling 
of W predicted by Eq. (J2J) with the exact values of 6 and y derived in Sec HU In Fig. |U we 
show the Log-Log plot of Wf(p) as function of t g(p) in the two limiting p values for fixed 
L. The results are in agreement with our exact results (See Sec |HJ and our scaling ansatz 
(Eq .(HZD). 

V. CONCLUSIONS 

We derive analitically the p dependance in the scaling behavior in two CGMs named 
BD/RD and EW/RD. Exact scaling exponents are derived and are in agreement with pre- 
viously conjectured values . To our knowledge these exact scaling behaviors were not 
analytically derived before. This derivation allows us to compute the scaling behaviors of 
the coefficients of the continous equations that describe their universality classes. We in- 
troduce the CEW/RD model to show that not all CGMs can be represented by an unique 
universality class. The CEW/RD is an EW in the limit p — > 1 while in the other limit p — > 
it followings the strong coupling behavior the KPZ equation. Our simulation results are in 
excellent agreement with the analytic predictions. 
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